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INTRODUCTION 

Quantifying entanglement in multipartite systems is 
a fundamental issue in Quantum Information Theory. 
However, for systems with more than two parts, very 
little is know in this respect. A few useful entanglement 
measures for pure states of 3 or 4 qubits have been inves- 
tigated 0, y| , but one is still far from a complete un- 
derstanding. Furthermore, for system of up to 4 qubits, 
a complete classification of entanglement patterns and of 
corresponding invariants under local filtering operations 
(also known as SLOCC, Stochastic Local Operations as- 
sisted by Classical Communication) is know jj, |5j . Kly- 
achko [a, proposed to associate entanglement (of pure 
states) in a /c-partite system (or perhaps , one should 
say "pure fc-partite" entanglement) with the mathemati- 
cal notion of semi-stability, borrowed from geometric in- 
variant theory, which means that at least one SLOCC 
invariant is non zero. For such states, the absolute val- 
ues of these invariants provide some kind of entanglement 
measure. However, even for system of k qubits, the com- 
plexity of these invariants grows very rapidly with the 
number of parts. For k = 2, they are given by simple lin- 
ear algebra The case k = 3 is already nontrivial but 
appears in the physics literature in [T^j an d boils down 
to a mathematical result which was known by 1880 [i"lj |. 
The case k = 4 is quite recent , and to the best of our 
knowledge, nothing was known for 5-qubit svstems|22j. 

Our main result is a closed expression of the Hilbert 
series of the algebra of SLOCC invariants of pure 5-qubit 
states. This result, which determines the number of lin- 
early independent homogeneous invariants in any degree, 
was obtained through intensive symbolic computations 
relying on a very recent algorithm for multivariate residue 
calculations. We point out a few properties which can be 
read off from the series, and determine the simplest in- 
variants, which are of degree 4 and 6 in the component 
of the states. 



HILBERT SERIES 



is H = 1/® 5 , which will be regarded as the natural repre- 
sentation of the group of invertible local filtering opera- 
tions, also known as reversible stochastic local quantum 
operations assisted by classical communication 

G = Gslocc = SL(2,C) x5 , 

that is, the group of 5-tuples of complex unimodular 2x2 
matrices. We will denote by 
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a state of the system. An element g = ( k g{) of G maps 
|*) to the state 

I*') =g|*) 
whose components are given by 

A' — V* 1 „Jl2 J23 734 745 75 A (-1 \ 
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We are interested in the dimension of the space Id of all 
G-invariant homogeneous polynomials of degree d = 2m 
(I4 = for odd d) in the 32 variables A^i^m^. 

It is known that it is equal to the multiplicity of the 
trivial character of the symmetric group ©2m m the fifth 
power of its irreducible character labeled by the partition 
[m, m] 

dimZ d = (x 2m |(x mm ) 5 ). (2) 
The generating function of these numbers 
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is called the Hilbert series of the algebra I = Q) d Id- 
Standard manipulations with symmetric functions allow 
to express it as a multidimensional residue: 
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2ttiui J 2niU5 B(u; t) 
where the contours are small circles around the origin 
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Denote by V = C the local Hilbert space of a two 
state particle. The state space of a five particule system 
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TABLE I: Coefficients of P(t) 



12 (Li, . . . , L$). These 17 polynomials are called the pri- 
mary invariants. 

The numerator should then describe the secondary in- 
variants, that is, a set of 3014400 homogeneous polyno- 
mials (1 of degree , 16 of degree 8, 9 of degree 10, 82 of 
degree 12 etc) such that any invariant polynomial can be 
uniquely expressed as a linear combination of secondary 
invariants, the coefficients being themselves polynomials 
in the primary invariants. 

This picture, which is the simplest kind of description 
to be expected, is far too complex for physical applica- 
tions. The best that can be done is to use the Hilbert 
series as a guide for finding explicitly a small set of rea- 
sonably simple invariants, in particular, the primary in- 
variants of lowest degrees. We have computed the first 
primary invariants, those of degree 4 and 6, using meth- 
ods from Classical Invariant Theory (see below). 



and 



THE SIMPLEST INVARIANTS 



B(u;t) = J][ (1 
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Such multidimensional residues are notoriously difficult 
to evaluate. After trying various approaches, we even- 
tually succeded by means of a recent algorithm due to 
Guoce Xin [l^, in a Maple implementation. The result 
can be cast in the form 
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where P(t) is an even polynomial of degree 104 with non 
negative integer coefficients a n 
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P(t) = a 2k t 2k 



k=0 

given in table [IJ and 

Q(t) = (1 - < 4 ) 5 (1 - t 6 )(l - i 8 ) 5 (l - i 10 )(l - t 12 f. 

On this expression, it is clear that a complete descrip- 
tion of the algebra of G-invariant polynomials by genera- 
tors and relations is out of reach of any computer system. 
Nevertheless, inspection of the Hilbert series suggests the 
following kind of structure for this algebra. We know, 
since dim7i — dimG = 2 5 — 3x5 = 1 7, that there must 
exist a set of 17 algebraically independent invariants. The 
denominator of the series, which is precisely a product of 
17 factors, makes it plausible that these invariants can be 
choosen as five polynomials of degree 4 (to be denoted 
by D x D y , D z , D t , D u ), one polynomial of degree 6 (F), 
five polynomials of degree 8 (Hi, H2 ■ . ■ ,-#5), one poly- 
nomial of degree 10 (J) and five polynomials of degree 



Transvectants and Cayley's Omega process 

In order to apply the formalism of Classical Invariant 
Theory, a state |^) will be interpreted as a quintilinear 
form on C 2 (called the ground form) 



il j 2 2 ,*3 ,*4*5— 
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A covariant of / is a G-invariant polynomial in the coef- 
ficients ^2122232425 and the variables Xi, yi, Zi, U and itj. A 
complete set of covariants can be in principle computed 
from the ground form by means of the so-called Omega 
process (see for notations). Cayley's Omega pro- 

cess consists in applying iteratively differential operators 
called transvections and defined by 

(P,Q)«i-«« = tx^---^P(x',...,u')Q(x",...,u") 
where 



fL = det 
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Degree 4 

Regarding parameter, write / as a quadrilinear 

binary form in the variables yi, Zi, U and Ui 



f ~ ^-j] 232324^1 Z »2^3 M » 



It is known that such a quadrilinear form admits an in- 
variant of degree 2 (called Cayley's hyperdeterminant 
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US EH E3 ) which is a quadratic binary form b x in the 
variables x = {xi,x%). Hence, taking the discriminant 
of b x one obtains an invariant D x of degree 4. We re- 
peat this operation for the other binary variables and we 
obtain four other invariants D yi D z , D t and D u . Evalu- 
ating the appropriate Jacobians with a computer algebra 
system gives the algebraic independance of the five in- 
variants. 



Degree 6 

We obtain the primary invariant of degree 6 by a suc- 
cession of transvections. First, we compute a triquadratic 
covariant of degree 2 

-D22020 — \J,J) 

This covariant allows to construct a cubico-quadrilincar 
covariant of degree 3 

L-31111 — 1^22020, J ) 

which gives a triquadratic polynomial of degree 4 

-D22200 = (C31111; Z) 1 11 • 

Hence, one obtains a quintilincar covariant of degree 5 

£■11111 = (-D22200: /) 11100 . 

Finally, we find the invariant of degree 6 

F=(E llllu f) llni . 

By computing the Jacobian, one finds that F is alge- 
braically independent of D x , . . . , D u . 

CONCLUSION 

From the Hilbert series, it appears that the algebra of 
polynomial invariants of a five qubit system has a very 
high complexity. Furthermore, as is already the case with 
smaller systems 0,0,01; the knowledge of the invariants 
is not sufficient to classify entanglement patterns. In the 
case of four qubits or three qutrits, this classification can 
be achived due to hidden symmetries which have their 
roots in very subtle aspects of the theory of semi-simple 
Lie algebras (Vinberg's theory 0]). However, such sym- 
metries are absent in the case of 5 qubits. Then, the only 
known general approach for classifying orbits (entangle- 
ment patterns) requires the computation of the algebra 
of covariants, which is already almost intractable in the 
case of four qubits. It has 170 generators, which have 
been found [3 , but the description of their algebraic re- 
lations (syzygies) is definitely out of reach. However, a 
closer look at the 4-qubit system, reveals that the classi- 
fication of Verstraete et al 0, . can be reproduced by 
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TABLE II: Evaluation of SLOCC covariants for Osterloh and 
Siewert states (x means that the evaluation is not 0) 



means of only a small set of covariants. We hope that our 
results will allow the identification and the calculation of 
such a small set of invariants and covariants, sufficient to 
separate the physically relevant entanglement patterns, 
which are probably not so numerous. To illustrate this 
principle, let us consider a result of Osterloh and Siew- 
ert |20j . Having introduced a notion of filter which can 
be used to separate SLOCC orbits in the same way as 
covariants, these authors show that the four states 

= (111111) + 100000)) 
|$ 2 ) = 2 (I 11111 ) + 111100) + 100010) + 100001)) 

I $3) = -p fx/2 111111) + 1 11000) + 100100) + 100010) 
\/6 v 

+ 100001)) 

I $4) = (\/3 111111) + 1 10000) + 101000) + 100100) 
2\/2 V 

+ 100010) + 100001)) 



are in different orbits. As can be seen on Table [H] the 
orbits of these states are also distinguished by our covari- 
ants. 

Finally, the investigation of entanglement measures re- 
quires an understanding of invariants under local unitary 
transformations (LUT) [2l| . In a forthcoming paper, we 
will explain how to obtain LUT-invariants from SLOCC- 
covariants. 
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